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ABSTRACT 


The generation of gravity waves in a frontal zone is investigated 
with a numerical model. The solutions show that gravity waves can be 
generated in a frontal zone by the application of a heating function 
which varies sinusoidally with time, and that such waves propagate both 
ahead of and behind the frontal zone. The amplitude of the waves are 
directly proportional to the magnitude of the forcing function and to 
the period of sinusoidal oscillation. Even with large heating rates, no 


hydraulic jumps are formed. 
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I. INTRODUCTION 


Meteorologists involved in numerical weather prediction have for some 
time been filtering gravity waves out of their prediction models. This 
has been done because high frequency inertial-gravity oscillations could 
be generated which would give unrealistic pressure changes, obscuring the 
meteorological forecasts. Gravity waves may be meteorologically signifi- 
cant ‘however in explaining irregularities in the ionosphere, and in some 
mesoscale phenomenon in the troposphere. Hines (1960) has suggested that 
internal atmospheric gravity waves may account for many of the irregulari- 
ties observed in the lower ionosphere and that these waves may have their 
origin in the large energy regions of the lower troposphere. Gossard, 
Jensen and Richter (1971) have studied gravity waves of small amplitude 
with a high-resolution FM/CW radar sounder. Breeding (1972) investigated 
the creation of unstable regions in the atmosphere by internal gravity 
waves. Gossard (1962) has listed the following three mechanisms for 
generating internal waves in the troposphere: 

"1. Internal waves can be generated as standing waves in the lee of 
topographic features. 

"2. Gravity waves are often caused by a traveling boundary between 
two types of air of different density. This type of tropospheric wave 
train is short, lasting less than two hours. It is a fair-weather 
phenomenon, rarely associated with storm fronts. 

"3, Oscillations which last for many hours are generally associated 
with tropospheric storms and synoptic scale features. These oscillations 
of long duration are usually generated by storms and frontal disturt -nces 


but the presence of a stable layer in the lower atmosphere is requi.cd." 





Tepper (1952) postulated that gravity waves may be generated in the 
vicinity of a moving frontal zone if the zone moves into an area where 
a stable layer exists. As the waves move out they steepen into hydraulic 
jumps. He felt that large vertical motions caused by the jumps could 
lead to severe weather such as is found near squall ines Kurth (1971) 
studied hydraulic jumps in a rotating, continuously stratified atmosphere. 
Numerical simulation of squall line formation was carried out by Sasaki 
(1959) and Ogura and Charney (1962). One of the objectives of the present 
study was to determine if such effects could occur near a properly simulated 
frontal zone. 

Since the advent of satellite pictures and the technique of producing 
film loops showing "“weather-in-motion" it has become evident that high 
speed gravity waves move out from the face Srudcental zones (Fujita and 
Bohan, 1967). 

In this study, the non-geostrophic frontogenesis model of Williams 
(99722) is modified in an aiceane to generate gravity waves in a frontal 
zone. 

In section II, the Boussinesq equations are given and the modeling 
relations are introduced. In the model the time-dependent quantities 
are functions of y and z only. The basic wind deformation field is 
independent of time and height. The numerical procedures, boundary 
conditions and initial conditions are given in section III. Numerical 
results are presented and discussed in section IV. Conclusions and 


recommendations for further investigations are given in section V. 





IT... “BASIC EQUATIONS 


In this experiment the Boussinesq equations are employed and the 
domain is bounded by two rigid horizontal planes. The compressibility 
of the atmosphere, which is neglected in the Boussinesq Bae civeicion. 
should not be of qualitative importance since the density scale height 
in the atmosphere is much larger than the thickness of typical frontal 
zones. The replacement of the tropopause by a rigid horizontal surface 
will induce large errors in this region and prevents gravity waves from 
propagating out of the troposphere, but the resulting errors near the 
lower boundary should be small. Basic equation development follows the 
work of Williams (1972a). 

The hydrostatic Boussinesq equations can be written in the following 


form when the earth's rotation and horizontal eddy viscosity are included, 


_ z 

AV +I W)tEWV) EVO FE KXYV =A VV, OED 

| ge 4. v- (eV) + Ww) - A Je + G@ ) (Et?) 

VV + ge =O, Gums) 

22 = s0/e,, (11.4) 

where °, is a constant reference potential temperature; p, a constant 
reference pressure, k = R/(Cp) ; O= Tpo/oe 5 the departure of 

the potential temperature from O,, P= Cp © (p/p0) + ez, the pressure 


2 
function, as Tv, the horizontal diffusion of momentum, A,V 6, the 


horizontal diffusion of heat., and Q the heating function. 
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An exact steady solution to this set of equations is given by 


V=D (xi-yj) 

w= 0 

p= b= | -D? (x"4y79/ 2 | -f Dxy » (IT.5) 
O=0 


where D is a constant. Cartesian coordinates have been introduced 
and f is taken to be constant. The horizontal velocity is a field of 
pure deformation and this deformation is given by 2D. 

If departures from the fields (11.5) are independent of x, then they 
will remain independent of x. Thus, the dependent variables are subdivided 


as follows: 


| p x+ w(y.z,t) | ibe | dytv(y.2,t) | J 
= w(y,z,t) 
O52, 0) ; (II.6) 


= ® (x,y) pee a (Cy, 2 ae 


If these expressions are substituted into Eqs. (II.1)-(11.4), they 


Il 


g = < 
I 


become: 
) yu 
Be + Gee au + VE fer = An Ge a 
SZ + he), aw) , pe ae » (11.8) 
90 4 (ve) t3 x(wo)4 V 38 =A. 38 a. +Q 3: AGLE=O) 
= 510 
Bs Be =O : 
Se a Gra!) 
oz Oo 
Where V = -Dy. The boundary conditions are 
we=oO, z=0, H Pee (idtad) 


where H is the distance between the horizontal plates. 


If the vertical average of a Sa is defined as 


<( =H) Ae, 


TT 





and the hydrostatic equation (II.11) is integrated with respect to z 


and the vertical mean is removed, then 


7-1) Es 2 fede-(/ede)] suet) 


The y equation of motion may now be written in terms of TT — Gi : 
If the vertical average of (II.10) and the boundary conditions (11.12) 


are used, then 
e \v7 =o : 


This equation states that the total disturbance y mass flux is independent 
of y. If the other variables have proper symmetry it follows that <v > 


must be an’odd function of y which leads to 


{v) = O : (11.14) 
If the vertical average of (I1I.8) is taken, and (II.12) and (11.14) 


are used, then 


uu ' wie, paca. 
oY dy GOES) 


This result is subtracted from (11.8) which yields 


Sv tg (vv-<vv)) +& (wr) +2 E(x V) + 2  (n- <1?) 


oc 
+f (u- <u)) = An yy (II.16) 


Bas. (11.7), (11.9), (11.10), (11.13) and (11.16) form a complete set 


which can be solved by a pure marching process. 


ez 





III. NUMERICAL PROCEDURE 


The arrangement of variables and the finite difference approximations 
are the same as those used by Williams (1967). In order to close the 
problem computational boundaries must be introduced in y. Since the 
disturbance velocity should die out at a sufficient distance from the 
axis of dilatation, then 

aaa. ¥ > Zou = 0 oe Cisich e0)) 
However, there is appreciable inflow across these computational boundaries 
since V+ Y) = + DY. The quantities u and @ which are advected across 
the boundaries must be specified independent of the interior values if 
computational stability is to be maintained (Platzman, 1954). Thus, the 


following boundary conditions are used: 


Wl Gray 2,2 ,¢] ultiYbeays 2), 2, 0] 


; (Dele) 
O[+ (Y+Ay/2,z,t] = e[t(%+Ay/2), z, 6] 


The computational boundaries y = + Y are placed between grid points 
so that the above conditions are actually applied at y=+ (Yt+tAy/2). 
These conditions are replaced by time specified functions in some 
experiments. 

These boundary conditions developed a weak nonlinear instability 
near the boundaries but this was removed by introducing a forward step 
every /2 time steps. 


The initial temperature field, from Williams (1972a), is given by 


S) (y,2, o) = (36, /)2) (2-H) ma? (2 Ar) are tan (simhery ) CLs) 


de 





payee! 
Pe, 


‘= Tn i 

were ovata (50, 00, / se) - The quantity 06. Pe 
which is constant, is the initial stability and A is one-half the total 
horizontal temperature variation. The initial x component of the velocity 


is given by , 
u(y 2,0) = & ai e& HA) sech (og). (III.4) 


The initial ~field is obtained from the quasigeostrophic equation. 
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IV. NUMERICAL SOLUTIONS 


The finite-difference fous of Eqs. (11.7), (129), Gao) ees 
and (II.16) are eoneaen a pure marching process. i goqrinenge are 
conducted by varying boundary conditions and heating functions. The 
Williams model is allowed to run for four days by which time a frontal 
zone has stabilized as described by Williams (1972b). At four days, 
boundary conditions are changed or heating is introduced and the model 
runs for an additional four days. The parameters which describe the 
various experiments are given in Table I. 

The procedure for initializing the heating function @ is as follows: 
At the end of four days, the vertical velocity field is scanned at 
z=2.5 km (level three) to determine the area of maximum vertical velocity, 
which corresponds to the leading edge of the frontal zone. A region R 
is defined which is the region where the heating function is applied. 
This region varies in width and vertical extent depending on the 


experiment. The function is applied as follows: 


Q 
Q = 0 elsewhere. 


© s(t) Gty,z) am region R, 


Here C is the maximum heating added and g(t) takes the value 1.0, or 
the following sinusoidal variation 

g(t) = s in(ws 
where Wis the frequency. The function G(y,z) takes either the value 


1.0 for block heating or the value determined by the sine squared function 


Z arLy— y, (2) ] 


G(y,z) = s in width of heating zone 


4 


In later .experiments, a cooling function is applied to the upper lec 21s 


of the model to simulate evaporation. This cooling is only activai d 


ws 











-gu0z Te}UCAJZ oY WoAZ Aeme pooeTdsTq 
‘qUOZ TeVUOAZ DY UT SANdIQ = 





dL) 
IEEE 
W) 


Y uortser ut (7-AeP Po ee pe 

Ga? ee 000T) Vv ‘VY | aisle 

(spuodes JO S,Q0QT) deas ewtTL EEE 

(wy JO 8,00T+) YWIPTA TePoW EEE 

ce RO 
LSS 


risa 620 a0 


OL. 6 i) 18") 


ie lt 7] so] sin2] sinae| 361-3 31.53D 


re [2030 sin? | einere| ep 3-7 Lis ome Tt 


id 6 | 18 | 
iiss) 


6 | 18 
6 | 18 


| 


sinwt| 36} 3-7 v}13.12DI 


| 1} 6 120] 30| sin} sine 36! 3-7 Vi13.68D. 
| 116 4 20| 30} sin? | sinest|52| 3-7 |v 18.27D) 


1 20 30} sin 


6 sine t} 361 3-7 713.300] | | 


1a [Pe 


al 
is 


oe 


Ao a8 
| 6 | 18 | 


ee 
4 


(18 {| 1/6 {20| 30| sin¢} sinwt! 20) 3-7 “| 3.20F 


aos 
17 


sin¢ | sinwt| 36|8-20v) 7.69D 





zy 


19 
20 


sin | sinwt|36| 3-7 |v|28. 46D 
sinwt | 36| 3-7 


100 | sin2 | sinwt 36| 3-7 poner 


Variation of parameters 


TABLE I. 


16 





when there is heating and it is ere otherwise. The magnitude of the 
cooling is equal to 204 of the heating below. 

Run number one is a control run with no variation of parameters. 
Maximum vertical velocity and the approximate position of the frontal 
zone is shown in Fig. la. Departure potential Lane eaeene: which is 
B= O~ eG, where 38; J) 2 is the initial stability, is shown in 
Fig. lb. The divergent part of the wind at four y grid points is 
shown in Fig. 2a and Fig. 2b. Fig. 3a and Fig. 3b show the variation 
in the departure potential temperature in the control run and the same 
quantity in run number three where a stable layer is introduced at the 
beginning. In both cases, all waves are virtually damped out at the end 
of eight days. 

Run number two tests Tepper's hypothesis that gravity waves will be 
generated if a frontal zone moves into an area where a stable layer exists. 
A stable iisipes is introduced at the end of four days ahead of the frontal 
zone at the -Y boundary between z=1.5 km and z=2.5 km where do; Voz 
was set equal to 8°K ain ee No gravity waves or jump zones were noted. 

mis number three is the same as the control run except for the 
introduction of a stable layer at both the +Y boundaries at the beginning 
of the program. As in run number two, the layer is inserted between 
Z2=1.5 km and z=2.5 km and 38, /3z is set equal to 8°K km between 
these two levels. By the end of four days this stable layer had 
propagated throughout the model between levels z=1.5 km and z=2.5 kn. 

The presence of a stable layer throughout the model did not appear to 


generate gravity waves or create jump zones. 


7, 





run. no. 1 
after 8 days 


Z (km) — 
O-WNWA OOD O 





1800 1200 600 O -600 -1200 -1800 
i Y ( km) 


Fig. la. Vertical velocity w maxima and the position of the 
front (---). Run number 1. 


run no,1 
after 8 days 





1800 1200 600 0 
- Y (km) 


-600 -1200 -1800 


Fig. 1b. Departure potential temperature 6' at level z=1.5 Im. 
Run number 1. 
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Fig. 2a. Divergent part of the wind v as a function of time at four 
levels: z=.s Km (1), z=2.5 km (3), z=5.5 km (6) and 2=3.5 mmc 
Run number 1]. 
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Fig. 2b. Divergent part of the wind v as a function of time at four 
Meveless2—-5 km (1), 2=2.5 km (3), 2=5.5 km (6) and z=c.5 knee 
Run number 1. 
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Fig. 3a. Departure potential temperature ©' as a function of time a 
four levels: z=.5 km (1), z=2.5 km (3), z=5.5 km (6) and z=8.5 km (9) 
Runs number 1 and 3. 










run no. 1 





run no. 3 
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In run number four, a heating function of 10°K Ales is introduced in 
block form and is 120 km wide, with vertical extent from z=.5 km to 
z=4.5 km, and it does not vary with time. This roughly approximated the 
release of latent heat which would be typical of convective activity in 
a frontal zone. As can be seen in Fig. 4a and Fig. 4b an increase in 
vertical velocity and a definite wave pattern in the departure potential 
temperature 9' ahead of the frontal zone was noted. The gravity waves 
appeared to be damped as the distance from the frontal zone increased. 
An additional result of this experiment was that the heating in the warm 
air greatly strengthened the temperature gradient in the frontal zone. 
This coincides with the work of Eliassen (1959) concerning the formation 
of fronts in the atmosphere. 

In run number five the errr function took the shape of a sine 
squared curve with maximum heating occurring in:the center of a 300 km 
wide zone with vertical extent from z=.5 km to z=4.5 km. This gave a 
smoother heating function in the frontal zone. The results were similar 
to run number four except that the vertical velocities and gravity wave 
amplitudes were reduced slightly, as might be expected. 

In run number six the heating function was allowed to vary sinu- 
soidaly with time which approximates the build-up and decay of convcctive 
activity in the vicinity of a frontal zone. The frequency for a gravity 
wave with a wave length equal to the scale of the heating zone gave a 
period of 9.9 hours. As a result of this variation, vertical velocity 
Maxima appeared at Y#%+1300 km which can be seen in Fig. 5a. Other 
bands of vertical velocity relative maxima appear but only the maximum 
values are shown in this and subsequent figures. Departure potential 


temperature fluctuations were also noted at Y#+1300 km as can be 2n 
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run. no. 4 
after 8 days 


Z (km) —-> 
OxwmN Wh UO 10 O 





1800 1200 600 O -600 -1200 -1800 
t= Y ( km) 


Fig. 4a. Vertical velocity w maxima and the position of the 
front (---). Run number 4. 
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run no. 4 
after 8 days 





-600 -1200 -1800 


1800 1200 600 0 
. + T (km) 


Fig. 4b. Departure potential temperature 0 ' at level 
z=1.5 km. Run number 4, 
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in Fig. 5b. More significantly, the gravity waves generated in the 
frontal zone propagated in both directions from the front and extend to 
the boundaries of the model. This can be seen in Fib. 6 which shows 
definite wave propagation in the departure potential temperature field 
at various locations on the Y-axis. 

In run number seven the heating function was reduced to 376 cae 
in order to determine if the patterns found in run number six would still 
be evident. Gravity wave generation and propagation still occurred 
but Ane magnitude of the waves was greatly reduced. Vertical velocity 
values also decreased and the maximum values returned to the vicinity 
of the frontal zone. This also was the first run where the heating 
function was raised and extended from z=2.5 km to z=6.5 km. This slight 
modification was made to account for the fact that convective activity 
in cumulus clouds does not extend to the surface. 

In run number eight, the heating function was increased to 30°K devia 
in order to increase the nonlinear effects in the gravity waves. Vertical 
velocity maxima remained at Y2%+1300 km but the magnitude of the vertical 
velocities were greatly increased. Also the first evidence of the 
formation of standing waves was observed. Well defined bands of vertical 
motion maxima were also noted. 

The heating function was increased to 50°K ae in run number nine 
and as can be seen in Fig. 7a and Fig. 7b, regions of strong vertical 
motion -were noted throughout the model as well as corresponding variations 
in potential temperature. Standing long waves appear to form with a 


wave length of approximately 360 kn. 
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Fig. 5a. Vertical velocity w maxima and the position of the 
front (---). Run number 6. 
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Fig. 5b. Departure potential temperature @' at level 
z=1.5 km. Run number 6. 
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Fig. 7a. Vertical velocity w maxima and the position of the 
front (---). Run number 9. 
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Fig. /b. Departure potential temperature O' at level z=1.5 kn. 
Run number 9, 
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In order to determine what effect a reduction of the horizontal 
eddy viscosity would have on gravity waves, in run number ten Am and 
Ag were reduced to 7,000m* sec~ from 20,000 m2 sec-1, This reduction 
appeared to increase vertical velocities and led to standing waves of 
approximately 130 km, as can be seen in Fig. 8a and Fig. 8b. 

In runs eleven through thirteen, the heating function was held at 
30°K day~1 with A, and Ay held at 20,000m2 secv! , and the stable layer 
was mieered: In run number eleven the stable layer was removed entirely 
and the areas of maximum vertical velocities shifted from Y*+ 1300 km 
to Y£+ 1200 km. In run number twelve a strong stable layer with 

dO, Jd? equal to 12°K km71 was inserted between z=1.5 km and 

z=2.5 km. The maximum vertical velocities moved almost to the extremities 
of the model, away from the frontal zone. Run number thirteeen was a 
repeat of run number two where Tepper's assumptions were tested, but with 
the inclusion of the heating function. Fig. 9a and Fig. 9b show much 

the same results as previous experiments with the exception that behind 
the frontal zone, where no stable layer exists, the maximum w was moved 
closer to the frontal zone. An interesting result of run number thirteen 
can be seen in Fig. 10 at Y=-840km and Y=-1140 km where a long wave 
developed at z=8.5 km (level nine). This agrees with Alberty and 
Van Sickle (1969) who found evidence of a mesoscale wave between a 
frontal zone and a severe storm. 

Run number fourteen was another control run which includes all the 
changes made to the original model. It offers an opportunity to make 
detailed comparisons with eue number one. Table I contains details of 


parameter changes made between run number one and run number fourteen. 
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Fig. 8a. Vertical velocity w maxima and the position of the 
front (---). Run number 10. 
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Fig. 8b. Departure potential temperature ©@' at level z=1.5 km. 
Run number 10. 
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Fig. 9a. Vertical velocity w maxima and the position of the 
front (---). Run number 13. 
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Run number 13. 
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Fig. 10. Departure potential temperature @' as a function of tin it 
four levels z=.5 km (1), z=2.5 km (3), 2z=5.5 km (6) and z2z=6.5 cio 
Run number 13, | 
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Fig. lla and Fig. 11b show the now common profile of w and ©' distribu- 
tions. Fig. 12a and Fig. 12b show the divergent part of the wind v as 

a function of time for four positions on the Y-axis ahead of the frontal 
zone. A phase velocity was obtained from the propagation of the phase 
lines of the velocity field. In this case the phase speed was approxi- 
mately 25msec"l, This can be compared to Alberty and Van Sickles' 
mesoscale wave velocity of 14m sec-l from’ observed data. Fig. 13 shows 
an interesting phenomenon which was present in all the cases where 
gravity waves propagate away from the frontal zone. As the waves at 
level three (2.5 km) and level six (5.5 km) propagate, the ©' traces st 
these two levels diverge. At the point in time where they first begin 
tomdaverge, the @'' ee at level nine (8.5 km) decreases while the €” 
trace at level one (0.5 km) eee slightly. When the ©@' traces at 
level three and level six begin to converge again, a long wave appears 
to be generated at level nine, with a marked increase in ©'. The ©' 
traces at Y=-1140 km seem to suggest that this divergent/convergent 
pattern repeats itself. No conclusions will be drawn concerning this 
phenomenon at this time, but it should be an interesting area for future 
theory and study. 

In runs number 15 and 16 the frequency with which the heating 
function varies with time is altered. In run number 15 a frequency of 
277/52,000 is used corresponding to a maximum every 14.5 hours. This 
led to an unrealistic slope in the frontal zone and an outward displace- 
ment of maximum w, as can be seen in Fig. l4a. The departure potential 
temepratures exhibited wide fluctuations as shown in Fig. 14b and 
Fig. 15. In run number 16 a frequency of 2417/20,000 is used which ‘ed 


CO a Wave period of 5.5 hours. The results of this experiment wees 
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Fig. lla. Vertical velocity w maxima and the position of the 
. front (---). Run number 14. 
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Fig. llb. Departure potential temperature @' at level z=1.5 km. 
Run number 14. | 
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Fig. 12a. Divergent part of the wind v as a function of time at four 
Hevels: z=0.5 km (1), z=2.5 km (3), 2=5.5 km @) and z=8.5 knoe 
Run number 14. 
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Fig. 12b. Divergent part of the wind v as a function of time at four 


Mmevels: z=0.5 km (1), z=2.5 lm (3), z=5.5 km (6), and z=8.5 km (oe 
Run number 14, 
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Fig. 13. Departure potential temperature ©@' as a function orse- 
at four levels: z=0.5 km (1), z=2.5 km (3), z=5.5 km (6) and Z-6.50er 
(9). Run number 14. 
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Fig. 14a. Vertical velocity w maxima and the positicn of the 
front (---). Run number 15. 
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Fig. 14b. Departure potential temperature @ ' at level z=1.5 km. 
Run number 15. 
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Fig. 15. Departure potential temperature @' as a function of time it 
four levels: z=0.5 km (1), z=2.5 km (3), z=5.5 km (6) and 2-6, 5.9 
Run number 15. 
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quite similar to the results of run number seven in that gravity waves 


L the waves 


were generated and did propagate but, even with C=30°K day 
generated were of very small amplitude. The maximum vertical motion was 
located in the frontal zone and the departure potential temperature 
profile was similar to that of run number one. 

| It was expected that if the nonlinear effects were important, that 
the scale of the waves would become smaller as the waves moved out from 
the source region. In order to further examine these effects, the 
boundaries were moved from 1800 km to 2700 km in runs 17--20, These 
experiments in fact do show a decrease in scale as the waves propagate 
out. If this decrease comes from nonlinear effects thenit should be 
more pronounced in cases with larger heating. If run 19 with C=3°K 
day~1 is compared to run 20 with C=100° day~1, it can be seen from 
Table I that the maximum w for run 20 is less than 33 times the maximum 
w for run 19. If a nonlinear scale reduction occurred in run 20 then 
the maximum w in that experiment would be larger than the maximum w in 
run 19, times the ratio of the heating rates. This seems to indicate 
that the observed scale may not be caused by the nonlinear effects. 
Rather the observed scale reduction probably comes from the contraction 
in the basic deformation field. 

A comparison was also made with the results obtained by Kurth (1971) 
who found that in a rotating, continuously stratified atmosphere, 
hydraulic jumps will form if a certain relationship exists between the 
Froude number and re eee number. The Froude and Rossby numbers 
were determined for runs seven, nine, fourteen, and twenty, using the 
Froude number F=V/[D(g Oa Onde 1/2 ], and the Rossby number R, = V/fL, 


With the result that all four cases fell on the jump side of the 


40 





iit) cal curve F=AR,“, where A is a constant in the range of 6.5-/.5. 
Even though jumps should have occurred, it is evident from the previous 
discussion that they did not occur. The conclusion must be that gravity 
waves, by themselves, do Bs eee hydraulic jumps under the conditions 


Giethis study. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The generation of gravity waves in a frontal zone has been investi- 
gated numerically. It was determined that the presence of a stable 
layer ahead of a frontal zone was not sufficient to zenerate gravity 
waves or hydraulic jumps. It was found that the addition of a heating 
function on the face of the front strengthens the frontal zone and 
generates gravity waves, but that these waves do not propagate away 
from the frontal zone. It was found that if the heating function 
applied to the frontal face is allowed to vary sinusoidaly with time, 
that gravity waves will be generated at the front and will propagate 
in both directions from the frontal zone and induce vertical motion 
maxima both ahead of and behind the front. It is believed that such 
motions could cause convective Pcie and condensation on the less 
stable warm air side. The amplitude of the gravity waves and the 
magnitude of the vertical motion maxima varied directly with the heating 
function and ane period of oscillation. A stable layer need not be 
present for gravity waves to generate and propagate, although the 
existence and strength of the stable layer appears to affect the 
distance the vertical velocity maxima are displaced from the front. 

The vertical velocity maxima Palate vee with the scale of the model 
but these fluctuations appear to be caused by the contraction in the 
basic deformation field. The amplitude of the gravity waves generated 
did not appear to be large enough to induce nonlinear effects which 
might lead to the formation of hydraulic jumps or squall lines. Squall 


lines do not appear to be formed by a simple Tepper mechanism but ~ iy 


still be caused by topographic effects. 
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In future studies using Williams model, a smaller grid scale should 
be used in order to obtain the fine detail which may be missed using 
a course grid. The model would be improved if moisture were parameter- 
ized or if it were possible to fix the forcing function based on 
observational data. The latter is not yet pcssible due to inability 
to measure the heating distribution. The model should be modified 
by replacing the rigid plate upper boundary with a boundary which would 
allow gravity waves to propagate out of the troposphere. It is also 
suggested that cooling be included both above and below the heating 
region. 

Possibilities for further study using observational data would be 
to determine if a correlation exists between stable layer strength and 
displacement of severe storms from a frontal zone, and also whether 
divergence and convergence of potential temperatures curves between 
two intermedieate levels in the troposphere leads to long wave formation 


near the tropopause. 
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